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Abstract 

Preliminary  to  the  stiidy  of  the  structure  of  an  optical  image 
from  the  standpoint  of  electromagnetic  theory,  an  integral  representation 
is  obtained, for  the  electromagnetic  field  in  the  image  space  of  an  optical 
system.  This  representation,  which  is  not  restricted  to  systems  of  low 
angular  aperture,  is  in  the  form  of  an  angular  spectrum  of  plane  waves, 
and  is  closely  related  to  that  introduced  by  Luneberg  (19UU)  as  a  vector 
generalization  of  well-known  formulas  of  Debye  (1909)  and  Picht  (192$). 
It  is  shown  that  the  representation  has  a  simple  physical  interpretation 
in  terms  of  a  modified  Huygens-Fresnel  principle  which  operates  with 
secondary  plane  waves  rather  then  xd.th  secondary  spherical  waves. 
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!•  Introduction 

Practically  all  investigations  which  relate  to  the  structure  of  the 
optical  image  have,  up  to  now,  been  almost  exclusively  based  on  geometrical 
optics  or  on  the  scalar  diffraction  theory.  Valuable  as  these  methods  are  in 
connection  with  the  design  of  optical  systems  and  the  evaluation  of  their 
performance  ,  they  suffer  from  two  limitations,  namely:  (1)  they  do  not  give 
indications  of  the  vectorial  features  of  the  image,  i.e.,  the  state  of  polari- 
zation of  the  li^t  in  the  image  region  and  the  direction  of  the  energy  flow 
(Poynting  vector) j  and  (2)  they  cannot  be  easily  applied  to  the  problem  of 
determining  the  structure  of  the  image  in  systems  of  wide  angular  aperture, 
such  as  certain  photographic  lenses  and  some  modern  astronomical  cameras. 
It  would  appear  that  more  refined  methods  are  needed  for  the  treatment  of 
such  problems.  Development  of  more  powerful  methods  also  seems desirable  be- 
cause of  probleuB which  arise  in  the  related  field  of  optics  of  microwaves. 
Microwave  systems  frequently  employ  numerical  apertures  which  are  considerably 
larger  than  those  encountered  in  systems  transmitting  light,  and  it  is  not 
always  clear  to  what  extent  the  usual  methods  can  be  used  with  confidence  [cf . 
Bachynski  and  Bekef 1,  195^ •   Moreover,  in  microwave  optics  the  knowledge  of 
the  state  of  polarization  of  the  image  field  is  of  direct  interest. 

In  the  present  investigation  an  integral  representation  is  obtained 
for  the  electromagnetic  field  in  the  image  space  of  an  optical  systera,which 
makes  it  possible  to  determine  the  co^^)lete  structxure  of  the  image  in  systems 
of  low  as  well  as  high  angular  aperture.  Since  the  specification  of  the  exact 
boundary  conditions  idiich  the  image  field  satisfies  present  formidable  diffi- 
culties, the  representation  is  an  approximate  one,  but  it  may  be  expected  to 
predict  the  field  id.th  a  high  degree  «f  accuracy  at  all  points  which  are  not 
too  close  to  the  exit  pupil.  (The  distance  from  the  aperture  must  be  large 
compared  to  the  wavelength.)  This  is  the  case  of  the  main  interest  in  con- 
nection with  practical  applications.  Our  representation  is  intimately  con- 
nected with  that  of  Luneberg  (19Uii)  and  may  be  regarded  as  a  generalization 
of  the  well-known  representation  introduced  for  the  scalar  case  by  Debye 
(1909)  and  extended  by  Picht  (1925)»  "^he  present  derivation  appears  to  be 

YoT  review  of  the  literatvire  see  Wolf  (1951)  •  A  survey  of  recent  research 
on  the  fovindations  of  some  of  the  methods  is  given  in  Wolf  (1955)  • 
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simpler  than  Luneberg's  and  shows  clearly  the  nature  of  the  appraximations 
that  are  involved  when  the  solution  is  applied  to  problems  of  practical  in- 
terest. It  is  also  shown  that  the  representation  has  a  sinqsle  physical  inter- 
pretation in  terms  of  a  modified  Hiiygens-Fresnel  principle  which  operates  with 
(vectorial)  secondary  plane  waves,  rather  than  with  secondary  spherical  waves. 

In  a  subsequent  paper  the  results  xd.ll  be  used  to  determine  the 
structtire  of  electroBiagnetic  images  in  the  focal  plane  of  aplanatic  systeras 
of  low  as  well  as  high  angular  aperture, 

2«  An  integral  representation  of  the  image  field 

Let  S  be  a  monochromatic  point  source  of  angular  frequency  co,  situ- 
ated in  the  object  space  of  a  rotationally  symmetrical  system  and  let 

(2.1)  E(P,t)  -  ^e(P)e-^*  ,   H(P,t)  =  ah(P)e"^* 

denote  the  electric  and  magnetic  vectors,  at  a  point  P  in  the  image  space,  at 
time  t.  The  vectors  e  and  h  which  are  to  be  determined  are  in  general  complex, 
and  the  symbol  J^  denotes  the  real  part. 

Let  s  denote  a  unit  vector  along  a  typical  geometrical  ray  in  the 
image  field,  i.e*,  a  ray  which  reaches  the  image  space  throu^  the  exit  pupil 
of  the  system.  The  positive  direction  of  s  is  taken  in  the  direction  of  the 
propagation  of  the  light  (Figure  l).   We  shall  show  that  within  an  accuracy 

adequate  for  practical  purposes,  e(P)  and  h(P)  for  any  point  P  in  the  image 
space  which  is  not  too  close  to  the  plane  of  the  exit  pupil  may  be  expressed 
in  the  form 

(2.2)  e(P)  -  (  ( tit)  e^^'^da  ,     h(P)  -  f  [  ^e^^^'^  dA   . 

Here  r  is  the  position  vector  of  P,  A(s)  and  B(s)  are  certain  (generally  complex) 
vector  functions  of  s  which  are  specified  by  the  geometric  optics  field,  and  the 
integration  is  taken  over  the  solid  angle  -TL  formed  by  all  the  geometrical  rays 
which  reach  the  image  space.  The  symbol  k  is  the  wave  number  in  the  image,  i.e., 

/„  -N         ,      CO    2n      2n      , 

(2.3)  k.n--^-n^-nk^, 

o 

vhere   n  is  the  refractive  index  of  the  image  space,  c  is  the  vacuum  wavelength, 

.  X  and  X  are  the  wavelengths  in  the  image  space  and  in  the  vacuum  respectively, 

and  k  is  the  vacuum  wave  number. 
o 
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Figure  1«     nitistrating  the  notation 
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The  integrals  (2.2)  represent  the  image  field  in  the  form  of  an 
angular  spectrum  of  plane  waves*   A  representation  of  this  type  is  always 
possible  provided  that  the  unit  vector  s  spans  all  directions  including 
complex  ones  (cf.  Bouwkan?)  (195U),  U1-U2]*  These  later  directions  represent 
surfcice  waves  or  evanescent  waves.  For  our  purposes  it  will  be  sufficient  to 
include  only  the  real  directions  which  are  associated  with  rays  in  the  image 
field  of  the  instniment  (s  belonging  to-TI), 

Let  X,  y,  z  be  the  components  of  the  position  vector  r  (coordinates 

of  P)  with  respect  to  Cartesian  rectangular  axes  with  origin  at  the  Gaussian 

image  point  S,  of  S  and  with  the  positive  z-direction  at  right  angles  to  the 

plane  of  the  exit  pupil  and  pointing  away  from  it  (Figure  1).  Further,  let 

s  ,  s  ,  s  be  the  con^jonents  of  s.  Since  s  is  a  unit  vector,  aiKi  the  light 

is  propagated  into  the  image  space  (s_  >  O), 

z 


(2.1;)  8^  -   I  /l  -  s2-  s2 


The  element  dfi-  of  the  solid  angle  is  given  by 
(2.5)  dH     . 


ds  ds 
3C  y 


so  that  the  first  integral  in  (2.2)  may  be  written  as 


r 


(2,6) 


e(P) 


A 


A( s^,Sy)  exp  [ik( s^x+Syy+s^jZ)! 


ds  ds 
«z 


there  is,  of  course,  a  strictly  similar  expression  for  h(P)« 
This  follows,  on  transforming  the  usual  expression 


(I) 


dfl  -  sin  9  cos  «  d/ 


where  (6,(!0  are  the  spherical  polar  angles.  We  have  s  =  sin  0  cos  0, 

s  «  sin  ©  sin  {^,  s  =  cos  6.  Hence 
y         ^*     z 

da         ds 


(II) 


ds^  ds   - 

X  y 


X      X 

ds    3s 


ded(2f 


sin  ©  cos  ©  ded^  =  s  sin  ©  d©dj^  , 


and  (2.5)  follows  from  (I)  and  (ll)< 
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To  determine  the  vectors  A  and  B,  consider  first  the  behavior  of  the 
integral  (2#6)  when  the  distance 


(2.7)         r  -  |?|   -   l/r^-  -2-  -^ 


X  +  7  +  2 


of  the  point  of  observation  P  from  the  origin  of  coordinates  S,  is  very  many 
wavelengths.  Let  u(tt  ,  u  ,  u  )  be  the  unit  vectors  in  the  direction  S-P,  i.e.. 


(2-8)     va^ .  5 ,  u  «  ^ ,  u^  -  ^  -  yr: 


"y   r  '    a   r     f~   ^"  y 


Using  (2*7)  and  (2«5),  the  integral  (2*6)  may  be  written  in  the  form 


(2.9)  e(P)     - 

where  fl 


-g £-  Mqpfikr  ^(s^,Sy}  ^3c*^y)]  ds^dSy.    , 


(2.10)         i^(\>SyJ  \*\)     =    Vx*  VV  ^z^z 

.VVV  IV^-VV  lr"V\     • 


^x 


Now  when  kr  »  2nrA  »  1>  th«  exponential  term  in  (2.9)  varies  rapidly  as  (s  ,3  ) 

X  y 

explores  the  domain  of  integration, and  the  real  and  imaginary  parts  wOl  change 
sign  many  times.  In  consequence  the  contribution  to  the  integral  from  the  dif- 
ferent elements  of  the  domain  of  integration  virtually  cancel  outj  the  only 

appreciable  contribution  comes  from  the  elements  (if  any)  centered  on  (s  ,s  ) 

X  y 

in  the  domain  of  integration  for  which  0  is  stationary.  This  is  the  essence 
of  the  principle  of  stationary  phase  which  may  be  used  to  find  the  asymptotic 
approximation  to  (2.9)  for  large  values  of  kr.  This  approximation  is  derived 
in  the  Appendix  at  the  end  of  this  paper.  It  is  shown  there,  that  when  z  <  0, 
i.e.,  when  P  lies  on  the  same  side  of  the  focal  plane  as  the  exit  pupil,  then 

(2.11)     e(P)  -  -  g  A(-u^,-Uy)  ^  +  &i^r       when  -  u  C  A  , 

-  ^(^)^^  when  -u<a   . 

The  symbols  £  and  ^  stand  for  ^belongs  to'  and  »does  not  belong  to*  respectively, 
i.e.,  the  first  line  in  (2.11)  applies  when  the  direction  PS^^  is  parallel  to  one 
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of  the  geometrical  rays  which  reaches  the  image  space,  and  the  second  line  applies 
when  there  is  no  such  ray.  When  z  >  0,  i.e.,  when  P  lies  on  the  side  of  the 
focal  plane  opposite  to  that  of  the  exit  pupil,  one  has,  in  place  of  (2.11), 

(2.12)  e(P)  -  2|  J(^^y  e___  ^  g^(  1.)  /        ^^^^    -  ^^   ^ 

.  di^)^^^  when  u^a   . 

Now  the  exact  boundary  conditions  which  the  field  vectors  satisfy 
cannot  be  determined  in  practice.  However,  if  the  linear  dimensions  of  the 
aperture  are  large  compared  to  the  wavelength,  and  if  P  is  not  too  close  to 
the  apertttre,  a  good  approximation  to  e  and  h  is  obtained  by  inposing 
the  approximate  Kirchhoff's  boundary  conditions,  i.e.,  the  conditions 

(2.13)  e(p')  -  ^^^P'),   h(p')  =  h^^^P*)     vhen  T*cJL      , 
t(P*)  "  0  ,       h(p')  -  0  when  pV '^   » 

where  e^  (P  ),  h   (P  )  are  the  values  of  the  xinperturbed  incident  field  at  a 
point  P  in  the  plane  of  the  exit  pupil  <M  ,     If  we  neglect  terms  of  order 
(kr)"^'^  we  have  from  (2.11)  and  (2.13) 

(2.1U)    t(3^,Sy)  -  -1^  t(i)(p')  r'  e^^'  , 

where  P  is  the  point  in  the  direction  u  =  -s  from  S,,and  r  »  X  is  the  distance 

t  t     -'■ 

S^P  .  Within  the  accuracy  here  in  question  P  may  evidently  be  also  identified 

with  the  pod-nt  in  the  aperture  through  which  the  ray  with  direction  of  propaga- 

tion  s  ■  -  u  passes. 

t 
Except  for  points  P  in  the  immediate  neighborhood  of  the  edge  of  the 

aperture,  e   (P  )  is,  to  a  good  approxiraation,  given  by  [cf.  Born  and  Wolf, 

1958,  Chapter  III] 

(2ol$)        e^^^(p')  -    ^     exprik^^(p')]   , 

In  a  typical  case  in  light  optics  r  >  10  cm,  X  ~  5.10  "^cm  so  that 
kr  »  2nr  A  ^  10  • 
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where  R,  and  R-  are  the  principal  radii  of  curvature  of  the  geometrical  wave 

1      2|      I  .^ 

front  through  P  ,  (S*  (P  )  is  the  eikonal  function  and  a  is  a  (generally  complex) 

vector  which  is  constant  along  the  ray  and  perpendicular  to  it.  We  shall  call 

a  the  strength  factor  of  the  ray  (s  ,s  ).  With  a  suitable  choice  of  the  phase 

constant  of  a,  (S  (P  )  represents  the  optical  path  from  the  source  point  S^  to 

p'»  From  (2.m)  and  (2.15) 

(2.16)  A(s^,3y)  -  -^    -^  exp[ik^(^(p')^3]    . 

We  may  express  (2.1U)  in  a  somewhat  different  form  which  involves  the 
wave  aberration  function  of  the  system.  Let  <rbe  the  Gaussian  reference  sphere, 
i.e.,  a  sphere  ^diich  has  its  center  at  the  Gaussian  ijnage  point  S^  of  S^  and 
which  passes  through  the  center  E  of  the  exit  pupil.  F\rrther,  let  W  be  the  geo- 
metrical wave  front  through  E.  If  F  and  G  are  the  points  in  which  P  S^  inter- 
sects (y   aixi  W  respectively,  we  have  (see  Figiire  2), 

(2.17)  Si?)   +nr'  -  [S^gJ  +  [gf]  +  |fsJ   , 

where  the  brackets  [••••3  denote  optical  lengtha.  Now  since  G  and  E  lie  on  the 
sarae  wave  front,  [s  g]   «  fs  eI  .  Also  IfsJ  =  nR  where  R  is  the  radius  of  the 
reference  sphere.  Hence 

(2.18)  <?(P*)  +  nr'   »    n(G  +  I  )     , 
where 


(2.19)  C     =    i  [s  e1  +  R     =    constant 

and 

(2.2C)  $     »    I  (s^,Sy)     -    i   [gf] 


is  the  aberration  function  of  the  system.  Also,  since  the  principal  centers 
of  curvature  of  the  geometrical  wave  front  lie  in  the  region  of  the  optical 
ijaage,  R-,  R^  and  r  can  only  differ  by  small  relative  amounts,  so  that  the 


The  aberration  is  usxially  measured  along  the  ray  through  F  rather  than  along 
the  radial  direction  ff ,  •  For  practical  puiposes  the  two  definitions  may  be 


taken  as  equivalent  [cf.  Wolf,  1952J. 
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Figure  2»     Definition  of  the  aberration  function. 
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term  r  /y^  y^    ^   (2.16)  may  be  replaced  by  unity.  If  we  also  substitute 
from  (2.18),  (2.16)  becomes 

(2.21)    A(s^,Sy)  -  -gt(s^,Sy)  exp[ik{c  +  |  (s^,Sy)}|  , 


* 


and  the  diffraction  integral  (2.6)  finally  takes  the  form' 

(2.22)    e(x,y,2)  >  -  g  e^^  I  I  — ^-JL.  explik  ^  J(s^,s  )+s^x  +  s  y  +  s^^zMds^ds  , 


1  1  ^  X    ^"^  1^^  ^^^  '^''^^  *^''''  *  ^^  ^  ^^'^]' 


111  a  strictly  similar  manner  the  corresponding  expression  for  h  is 


obtained: 

(2.23)    h(x,y,z) 


JJ  X  y 


Here  b  is  the  magnetic  strength  factor,  i.e.,  the  vector  which  bears  the  same 
relation  to  the  incident  magnetic  field  h^'''^  as  a  bears  to  e^   |cf.  (2.15)3 • 
It  follows  from  Maxwell's  equations  that  a  and  b  obey  the  relation  [cf.  Born  and 
Wolf  (1958),  <§  3.1  eq.  (12a)  and  (l5a)] 

(2»2lt)  b  »   /eTP"  s^a  , 

where  e  and  |i  denote  the  dielectric  constant  eind  the  magnetic  permeability  of 
the  image  space,  and  the  relation  n  =  yiJI  is  assumed.  In  our  system  of  \inits 
(Gaussian)  \l  differs  inappreciably  from  unity,  but  is  retained  here  for  the 
sake  of  symmetry. 

The  formulas  (2.22)  and  (2.23)  represent  the  solution  of  our  problem. 
They  express  the  image  field  as  a  superposition  of  plane  waves  of  different 
directions  of  propagation.  The 'conqslex  an5)litudes'  a,  b  of  each  constituent 
plane  wave  can  be  derived  from  a  ray  trace,  taking  into  account  the  state  of 
polarization  of  the  field  along  each  ray.  The  phase  function  J,  which  is  a 
measure  of  the  aberration  of  the  system,  may  likewise  be  determined  from  the 
ray  trsice.  It  follows,  in  view  of  (2.21;),  that  our  solution  satisfies  the 

Since  the  origin  of  the  phase  is  arbitrary,  we  may,  of  course,  set  C  =  0. 
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homogeneoTis  Maxwell's  equations 

(2.25)  curl  h  +  ik  e  e  =  0  , 

curl  e-ikjih  =  0  , 
o 

throughout  the  whole  image  space,  and  it  is  seen  from  (2»12)  that  it  does  obey 
the  vectorial  form  of  the  Sommerfeld  radiation  condition  (cf.  C.  ^ffille^  (19$7), 
p»  13-1^3  in  the  half -space  z  >  0.  ^nce  the  apprciximate  Kirchhoff  boundary 
conditions  were  used,  our  solution  cannot  be  expected  to  represent  the  field 
adequately  at  points  in  the  immediate  neighborhood  of  the  exit  pupil.  However, 
when  the  aperture  is  large  compared  to  the  wavelength,  and  when  the  point  of 
observation  P(x,y,z)  as  well  as  the  Gaussian  focus  S,  are  both  situated  at 
distances  from  the  aperture  that  are  also  large  conq^ared  to  the  wavelength, 
(2o22)  and  (2 •23)  may  be  expected  to  give  a  good  approximation  to  the  field. 

3,  A  ptysicetl  interpretation  of  the  diffraction  integrals 

The  integrals  (2,22)  and  (2 .23)  have  a  sinqsle  physical  interpretation 

which  is  clearly  brought  out  when  comparing  them  with  the  corresponding  integrals 

based  on  the  application  of  the  Hviygens-Fresnel  principle.  According  to  the 

Huygens-Fresnel  principle,  each  point  G  of  the  wave  front  W  may  be  regarded  as 

a  center  of  a  secondary  disturbance  which  is  propagated  as  a  spherical  wave  of 

-4i)     -*(i)  — — — ..^— — 

complex  amplitude  proportional  to  e^  "^(G),  h^  '(G).  The  total  field  at  P  is 

obtained  by  adding  the  effects  of  all  the  spherical  waves  which  reach  the  image 

space  through  the  exit  pupil.  Thus  the  Huygens-Fresnel  approximation  (suffix 

m)   for  the  field  is 

(3.1)         ^(P)  -    jj   ^^^(G)  ^  K(G)  dW  , 


^(^> 


^ . X     ikd 
h^^^(G)  ^-  K(G)  dW 


where  d  is  the  distance  GP  and  K  is  the  usual  inclination  factor,  which, for 
sriiall  obliquities,  has  the  constant  value  -ik/2n.  The  integration  extends  over 
the  part  of  the  wave  front  which  approximately  iills  the  exit  pupil  (Figure  3(a)). 
In  the  present  formulation,  the  field  is  given  by 
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(a)    Huygens'    construction   relating 


to   the   integrals    (3.1). 


Spherical     Secondary    Waves 


(b)    A  construction    relating  to 
the    integrals    (3.2), 


Plane    Secondary  Waves 


•j^,*^    iks-r        di,-»,    iks.r 
A(s)e         ,       B(s)e 


Figure  3.  Physical  interpretation  of  the  diffraction  integrals. 
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(3.2)  t(P)  -    )}      Ait)   e^*^  dH 

n 
h(P)  =       1 1    B(i)  e^^-^  dn 


I 

Jf  B(^ 


a 

i.e.,  as  superposition  of  plane  waves  of  complex  ariplitudes  A(s),  B(s),  propa- 
gated  in  the  s-directions  coinciding  with  the  directions  of  all  the  geometrical 
rays  that  reach  the  image  space.  To  each  plane  wave  there  belongs  a  plane  wave 
front  which  is  tangential  to  the  wave  front  W  at  the  appropriate  point  G  (see 
Figure  3(b)).  Thus  each  point  of  that  portion  of  the  wave  front  W  which  approxi- 
mately fills  the  exit  pupil  may  be  regarded  as  giving  rise  to  a  secondary  plane 
wave,  and  the  field  at  P  is  obtained  by  adding  the  effects  of  all  these  secondary 
waves.   If  terms  of  order  (1/kr)  '  are  neglected,  the  vector  aii5)lit\xdes  in  the 
two  constructions  are,  according  to  (2.11),  related  by 

(3.3)  e(i)(G)  .  -g  A(^)  e^   , 

-(i).^.      2n  -,^  e-^^ 
h'  '(G)  «  -  j^  B(s)  — jr—   , 

where  r  is  the  distance  S,Q. 

Our  diffjraction  integrals  (2.22)  and  (2o23)  are  closely  related  to  the 
diffraction  integrals  of  Luneberg  (19UU)  and  may  be  regarded  as  generalizations 
of  the  integral  introduced  by  Debye  (1909)  to  represent  a  diffracted  spherical 
scalar  wave.  A  generalization  of  Debye »s  representation  to  non-spherical  (scalar) 
waves  was  carried  out  by  Picht  (1925) j  this  case  was  also  treated  recently  by 
Focke  (1956),  who  employed  essentially  the  same  method  as  used  in  the  present 
paper. 

Luneberg' s  integrals  maybe  derived  from  our  solution  as  follows:  We 
have  from   (2.19)  and  (2*20), 

(3.1;)        k[c  *l]   -  k  /[S^E]  .  [gf]  *   [FSjj 

Let  N  be  the  foot  of  the  perpendicular  from  S,  onto  the  ray  through  G  (Figure  h)< 
In  the  asymptotic  approximation ,  (gfJ  +  ^^t1  ™^y  evidently  be  replaced  by  the 
optical  length  (gnJ  measured  along  the  ray.  Further,  since  G  and  E  lie  on  the 
same  wave  front,  JS  E]  »  [s  g]|,  and  we  may  write  in  place  of  (3.I4) 
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Figure  h;     Transition  to  Luneberg's  diffraction  integrals. 


(3.5) 


[C  +  1}       -     k^|^[s^G]   ^    [ON]] 


"     KM 


Let  MS   now  introduce  in  place  of  the  \anit  vector  s  the  'ray  vector' 


n  s  , 


and  let  p,  q  and  in  be  its  coinponents,  i.e., 
(3.7)  P 


ns^,   q 


nsy. 


m 


ns. 


I  y  n  -  p  -  q   I   • 


Then  JS  NJ  becomes  a  function  of  p  and  q,  and  of  the  position  of  the  source 
point  S  }  it  is  the  mixed  characteristic  function  W(S  j  p,q)  of  Hamilton.  The 
strength  factors  a  and  b,  regained  as  functions  of  p  and  q  rather  than  of  s^ 
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and  s  ,  vd.ll  be  denoted  by  a  and  b  respectively,  and  the  corresponding  domain 
of  integration,  formed  by  the  (p,q)-coDiponents  of  all  the  geometrical  rays  which 
reach  the  image  space,  will  be  denoted  hy  Sl  ,     The  integrals  (2,22)  and  (2o23) 
become 


(3*8)    e(x,y,z) 

i3>9)  h(x,y,z) 

In  place  of  (2»2li)  we  now  have 

(3.10)  b 


2l£i22  expfik:^[w(S^}p,q)   +  px+qy+mz  j  Idpdq   , 
^  l^'^     exprik^[w(S^jp,q)  +  px+qy+mzj    |dpdq  • 


1    ^ 


A? 


The  formulas  (3»8)  and  (3»9)  are  formally  identical  with  Luneberg's 
diffraction  integrals.   They  were  obtained  by  Luneberg  not  as  an  approximate 
solution  of  the  true  physical  problem,  but  rather  as  an  exact  solution  to  an 
idealized  problem,  namely  the  following:  to  find  a  solution  of  the  homogeneous 
Maxwell's  equations,  valid  throughout  the  whole  space,  which  has  a  prescribed 
behavior  at  infinity.  In  spite  of  the  mathematical  elegance  of  Luneberg's 
formulation,  we  prefer  to  regard  our  integrals  as  an  approximate  solution  of  the 
true  physical  problem,  valid  at  distances  from  the  aperture  that  are  large  com- 
pared to  the  wavelength.  Our  analysis  shows  that  the  integrals  (2 .22)  and 
(2,23)  [or  (3«8)  and  (3 •9)]  are  then  effectively  equivalent  to  the  solution 
obtained  by  a  direct  application  of  the  (vectorial)  HTjygens-Fresnel-Kirchhoff 
diffraction  theory.  However,  since  our  formulas  do  not  involve  an  inclination 
factor,  they  are  more  suitable  for  the  treatments  of  problems  involving  large 
obliqxiities,  e.g.,  the  problem  of  determining  the  stnicttire  of  an  optical  image 
in  a  system  with  a  high  angular  apertrire.  This  problem  will  be  disctissed  in 
another  paper. 
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Appendix 
Asymptotic  evaluation  of  the  integral  (2«9) 

We  derive  here  the  asymptotic  approximation  for  large  kr  to  the 

integral  (2 ,9),  viz., 

ff   Ms^,sJ 
(Al)  e     » 

where  XI 


s'^      exp[ikri2f(s^,Syj  ^>^y)]  ^^x^j     * 


(A2)  J^(^x'V  ""x'V     '    Vx*  VV  ""z^z       » 

(A3)  s^  =     I  >4-s^-  Sy  I     , 


(Ai*) 


■     yl-  u  -  u  • 


As  explained  in  Section  2,  the  asymptotic  approximation  may  be  obtained  by 
the  amplication  of  the  principle  of  stationary  phase.  The  principle  is  well 
known  in  connection  with  single  integrals;  its  generalization  to  double  inte- 
grals of  a  wide  class  of  practical  interest  was  carried  out  by  Focke  (1951i), 
who  determined  the  whole  asymptotic  expansion  of  these  integrals.  Other  ways 
of  deriving  the  asyn^Jtotic  expansions  ware  discussed  by  Jones  and  EU.ne(1956), 
and  Chako  (1958).  It  follows  from  these  investigations,  that  for  large  kr, 
(Al)  reduces  to 


t  t 


A(s^,sJ    r.  .,  I  I     n   «;,1^3/2 


Here  s  ,  s  are  the  values  of  s  and  s  for  which  ^  is  stationary  within  the 
x'  y  X     y 

domain  of  integration,  i.e.,  they  are  the  roots  of  the  equation 

(^«  i-  ■  i-  ■  °  » 


"^Strictly  speaking,  the  »aD?)litude  function'  in  the  integrals  considered  by  these 
authors  is  a  real  scalar,  not  a  conqjlex  vector  as  in  (Al),  but  the  generalization 
to  the  present  case  is  quite  trivial. 
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further , 
(A7) 


(A8) 


A  = 


2^   2 


da       ds 


L   when  A  >  0  , 

j  s  "^-1    when  A  >  0  , 

L   when  A  <  0  , 


>  0 


<  0 


and  the  prinre  in  (A?)  and  (A8)  denotes  values  at  the  stationary  point.  The 
summation  in  (A5)  is  taken  over  all  the  stationary  points. 

If  there  is  no  stationary  point  within  the  domain  X\  of  integration, 
then 


(A9) 


1  3/^ 


(AlO) 


Now  we  have  from  (A2)  and  (A3), 

s^                                   ds  3 

_£_-   a   U  +  U   ■^—   a   U  -  U   — 

ds       X   z  ds       X  Z  S 


with  a  similar  expression  for  d0/ds   .  It  follows,  if  (Aii)  is  also  used,  that 
0  is  stationary  when 


(All) 


u 


s  =  u  , 
X     x* 


t  I 

3   =   U  ,    S   =   U 

y    y    "    " 


u 


z 
I 


u_ 


if  u^  >  0,  (a) 


if  u„  <  0.  (b) 


Thus,  if  u  >  0,  there  is  one  stationary  point  which  contributes  to  the  asyn^^totic 

approximation  of  e,  provided  that  v.   lies  within  JTj  if  u  <  0,  the  same  is  true, 

provided  that  -u  lies  in.fl*     If  u  (or  -u  ,  as  the  case  may  be)  does  not  lie  in 

z      z  _^ 

jfi ,  ^  is  not  stationary  within  the  domain  of  integration,  and,  by  (A9)  e  is  then 
of  the  order  of  (kr)  '  . 
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From  (AlO)  and  the  corresponding  relation  for  30/3s  .  we  have,  on  using 
(All), 


where  the  upper  sign  is  taken  when  u„>  0  and  the  lower  sign,  when  u„<  0«  The 

z  z 

expressions  (A7)  and  (A8)  become 
(A13)  A  -  ^   , 

(Allj)  j  »  +  1       according  as  ti  <  0  • 

We  also  have 

(A15)  ^^Sx'4^  '    ^^-  V  -  "V^  ' 

(A16)  ^(^'V  ^x*V  '-■'■* 

the  upper  or  lower  slga  again  being  taken  according  as  u^  ^    0. 

From  {A$)  and  (A9)  we  thus  finally  obtain  the  following  asyi^totic 
approximations  to  e  for  large  kr: 

(A17)  e     -      ||A(n^,Uy)^    +    0(^)  when    Ujj>  0    and    ucH, 

m      ^i^)  vhen    u>0    and    rxf.-^f 


(a) 


(b) 
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Boeing  Airplane  Company 

Seattle  Division 

Seattle  Uj,  Washington 

Attnt  E.T.  Allen,   Library  Supervisor 

Convalr,   A  Division  of  General  Dynamics  Corp. 

Fort  Worth,  Texas 

Attnt   K.Q.  Brown,   Division  Research  Ub. 

Convalr,   A  Div.   of  Qeneral  Dynamics  Corp. 

San  Diego  12,   California 

Attni  Mrs.   Dora  B.  Burke,   Engr.  Lib. 

Cornell  Aero.   Laboratory,   Inc. 

U1j55  Genesee  Street 

Buffalo  21,  New  York 

Attn:  Elma  T.  Evans,   Librarian 

Daljno  actor  Company 

A  Division  of  Textron,   Inc. 

1515  Industrial  Way 

Belmont,   California 

Attnt  Mary  Ellen  Addems,   Tech.  Lib. 

Dome  and  Margolin,   Inc. 
30  Sylvester  Street 
Westbury,   Long  Island,   N.Y. 
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Douglas  Aircraft  Company,  Inc. 
P.O.  Box  200 
Long  Beach  1,   California 
Attn:  Engineering  Lib.    (C-250) 

Douglas  Aircraft  Company,    Inc. 
827  Laphan  Street 
El  Segundo,   California 
Attnt  Engineering  Library 

Douglas  Aircraft  Company,   Inc. 

3000  Ocean  Park  Boulevard 

Santa  Monica,   California 

Attnt  Eq.  Sec.  Rof.  Files,  Eq.  Eng. 

Douglas  Aircraft  Corporation 
2000  North  Memorial  Drive 
Tulsa,  Oklahoma 
Attnt  Engineering  Library,   D-250 

Electronic  Communication,   Inc. 

St.  Paul  Street  at  University  Parkway 

The  Marylanda-  Apartments 

Baltimore  18,  Maryland 

Attn:  Dr.  D.D.  King,   Vice-President 

Electronic  Defense  Laboratory 
P.O.  Box  205 

Mountain  View,   California 
Attn:  Library 

Electronics  Research,  Inc. 

2300  N.New  York  Avenue 

Evansvllle,   Indiana 

Attn:  B.H.  Baldrldge,  Engr.  Supervisor 

Emerson  and  Cunlng,  Inc. 
869  Washington  Street 
Canton,  Mass. 
Attn:  W.R.  Cuming 

The  Qnerson  Electric  Mfg.   Company 

8100  Florissant  Avenue 

St.   Louis  21,  Missouri 

Attn:  Mr.   E.  Breslin,  Librarian 


Fairchild  Aircraft 

Division  Fairchild  Eng. 
Hagerstown,  Maryland 
Attnt  Library 


Farnsworth  Electronics  Company 
Fort  Wayne  1,   Indiana 
Attn:     Technical  Library 

Foderal  Telecommunication  Laboratories 
500  Washington  Avenue 
Nutiey  10,  New  Jersey 
Attn:  Technical  Library 

Oabriel  Electronics 
Division  of  Gabriel  Company 
135  Crescent  Road 
Needham  Heights  9U,  Mass. 
Attn:  Mr.   Steven  Oalagan 

Qeneral  Electric  Company 
Ithaca,  New  York 
Attn:     J.B.  Travis, 

Advanced  Electronics  Center 

General  Electric  Company 
Electronics  Park 
Syracuse,  New  York 

Attnt  Bessie  Fletcher,   Documents  Library 
Bldg.  3-lli3A 

General  Precision  Laboratory,   Inc. 

63  Bedford  Road 

Pleasantville,  New  York 

Attn:  Mrs.  Mary  G.   Herbst,   Librarian 

Goodyear  Aircraft  Corp. 
1210  Massillon  Road 
Akron  15,   Ohio 
Attnt  Ubrary  Vi/\20  Plant  A 

Granger  Associates 

Electronic  Systems 

966  Commercial  Street 

Palo  Alto,  California 

Attn:   John  V.N.  Granger,   Pres. 


Orxiiman  Aircraft  Engineering  Corp. 
Bethpaga,  Long  Island,   N.Y. 
Attnt     Mrs.  A.M.  Gray,   Librarian 

Engineering  Library,   Plant  No.  5 

The  Hallicrafters  Co. 
IjliOl  West  5th  Avenue 
Chicago  2li,  Illinois 
ATTN:   L.  La  Giola,   Librarian 

Hoffman  Laboratories,    Inc. 
37UO  South  Grand  Avenue 
Los  Angeles  7,    California 
Attnt  Engineering  Library 

Hughes  Aircraft  Co. ,   Antenna  Dept. 
Microwave  Laboratory,   Bldg.   12,   Ra  2617 
Culver  City,   California 
Attnt  M.D.  Adcock 

Hughes  Aircraft  Co. 
Florence  and  Teale  Streets 
Culver  City,   California 
Attn:  Dr.  L.C.   Van  Atta,   Assoc.  Dlr. 
Research  Laboratories 


Hycon  Eastern,   Inc. 

75  Cambridge  Parkway 

Cambridge,  Mass. 

Attn:  Mrs.    Lois  Seulowltz,   Tech. 


Ub. 


and  Airplane  Corp. 


International  Business  Machines  Corp. 

Military  Products  Division 

590  Madison  Avenue 

New  York  22,   New  York 

Attn:  Mr.  C.F.  McElwain,   General  Manager 

IntRrnational  Business  Machines  Co. 
Military  Products  Div. 
Oswego,   New  York 
Atnn:  Mr.  D.I.  Marr,  Librarian 
Dept.  1:59 

International  Resistance  Ccrapany 
bOl  N.  Broad  Street 
Philadelphia  8,   Pa. 
Attn:  Research  Library 

Jansky  and  Bailey,  Inc. 
1339  Wisconsin  Avenue,   N.W. 
Washington  7,   D.C. 
Attn:  Mr.  Delmer  C.  Ports 

Dr.   Henry  Jaaik,   Consulting  Engr. 
298  Shames  Drive, 
Brush  Hollow  Industrial  Park 
Westbury,  New  York 

Martin  Katxin  and  Company,   Consultants 
711  llith  Street,  N.W. 
Washington  5,   D.C. 

Lockheed  Aircraft  Corp. 

Missile  Systems  Division 

7701  Woodley  Avenue 

Van  Nuys,    California 

Attn:   Van  Nuys  Lib.  Bldg.  902 

Lockheed  Aircraft  Corp. 
California  Division  Engr.  Lib. 
Dept.   72-25,   Plant  *-l.   Bldg.   63-1 
Burbank,    California 
Attn:  N.C.  Harnole 

The  Glenn  L.  Martin  Company 
Denver  Division 
Waterton,  Colorado 
Attn:  WilUam  V.   Foley, 

Test  and  Reliability 

Mail  No.   2-38 

The  Glenn  L.  Martin  Company 

Baltimore  3,  Maryland 

Attn:  Engr.  Ub.,  Antsnna  Design  Group 

Harylaixi  Electronic  Mfg.  Corp. 
5009  Calvert  Road 
College  Park,  Maryland 
Attnt  Mr.  H.  Warren  Cooper 

Mathematical  Reviews 
190  Hope  Street 
Providence  6,  Rhode  Island 
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The  W.L.  Maxson  Corporation 

li60  West  3Uth  Street 

Vew  York,  (J.Y. 

Attn:  Miss  Dorothy  Clark 

MoDonnell  Aircraft  Corp. 
Lambert  Saint-Louia  Municipal  Airport 
Box  ^6,   St.  Lovds  3,  Misaouri 
Attn:  R.D.  Detrlch,   Engr.  library 

McMillan  Laboratory,  InC' 

Brojmvlile  Avsnue 

Ipswich,  Maasachuaetts 

Attn:  Security  Officer,   Documant  Rm. 

Melpar,  InC 

3000  Arlington  Blvd. 

Falls  Church,   Va. 

Attni  Engineering  Tech.  Lib. 

Microwave  Developmant  Laboratory 

90  Broad  Street 

Babson  Park  57,  Mass. 

Attni  N.  Tucker,  General  Manager 

Microwave  Radiation  Company 
1?223  South  Hamilton  St. 
Sardena,   California 
Attn:  Mr.  Morris  J.   Ehrllch,  Pres. 

Naval  Industrial  Reserve  Aircraft  Plant 

P.O.  Box  5907 

Dallas,  Texas 

Attn:  Mr.  H.S.  White,  Ubrarian 

Horthrop  Aircraft,  Inc. 

Hawthorne,  California 

Attnt  Mr.  E.A.  Freitas,  Lib.  Dopt.  31Ji5 

North  American  Aviation,  Inc. 
12211:  Lakewood  Blvd. 
Downey,  California 
Attnt  Engr.  Lib.  U95-115 

North  American  Aviation,  Inc. 
Los  Angeles  International  Airport 
Los  Angeles  1:5,  California 
Attn:  Engineering  Technical  File 

Page  Communications  Engineers,  Inc. 
710  Fourteanth  Street,  Northwest 
Washington  5,  D.C. 
Attn:  Librarian 

Phileo  Corporation  Res.  Division 

Branch  Library 

U700  Wlssachlckon  Avenue 

Philadelphia  hk,   Pa. 

Attnt  Mrs.  Dorothy  S.  Collins 

Pickard  and  Burns,  Inc. 
21:0  Highland  Avenue 
Needhan  9U,  Mass. 
Attn:  Dr.  J.T.  deBettenoourt 

Polytechnic  Research  and  Dev.  Co.,  Inc. 
202  Tillary  Street 
Brooklyn  1,  New  York 
Attnt  Technical  Library 

Radiation  Engineering  Laboratory 
Main  Street 
Maynard,  Mass. 
Attn:  Dr.  John  Ruie 

Radiation,  Inc. 
P.O.  Drawer  37 
Melbourne,  Florida 
Attn:  Tech.  Library,  Mr.  M.L.  Cox, 
Assistant  Project  Engineer 

Radio  Corp.  of  America 

RCA  Laboratories 

Rocky  Point,  New  York 

Attn:  P.S.  Carter,  Lab.  Ubrary 

RCA  Laboratories 

David  Sarnoff  Research  Center 

Princeton,  New  Jersey 

Attnt  Ml8»  Fern  Cloak,  Lib.,  Res. 


(2jmo  Ramo-Wooldrldge  Corporation 
5730  Arbor  Vltae  Street 
Los  Angeles  U5,  California 
Attnt  Margaret  C.  Whltnah,  Chief  lib. 

Hoover  Microwaves 
9592  Baltimore  Avenue 
College  Park,  Maryland 

The  Rand  Corporation 

1700  Main  Street 

Santa  Monica,  California 

Attn:  Library  -  via:  SBAMA  Uaison  Office 

Rantec  Corporation 

Calabasas,  California 

Attn:  Qrace  Keener,  Office  Manager 

Raytheon  Manufacturing  Company 
Missile  Systems  Division 
Bedford,  Mass. 
Attn:  Mr.  Irving  Goldstein 

Raytheon  Manufacturing  Company 
Wayland  Laboratory 
Wayland,  Masa. 

Attnt  J.E.  Walsh,  Dept.  897U, 
Antenna  Design  Radar 

Raytheon  Manufacturing  Company 

Wayland  Laboratory 

Wayland,  Mass. 

Attn:  Miss  Alice  0.  Anderson,  Librarian 

Republic  Aviation  Corporation 
Farmlngdale,  Long  Island,  N.Y. 
Attn:  R.E.  Fidoten,  Engineering  Lib. 

Rhesm  Manufacturing  Company 
Research  and  Development 
9236  East  Hall  Road 
Downey,  California 
Attn:  J.C.  Joerger 

Ryan  Aeronautical  Company 
Lindbergh  Field 
San  Diego  12,  California 
Attn;  Library 

Sage  Laboratories,  Inc. 
30  Qulnan  Street 
Walthan  51i,  Mass. 

Sanders  Associates,  Inc. 
95  Canal  Street 
Nashua,  New  Hampshire 
Attn:  N.R.  Wild,  Library 

Sandia  Corporation,  Sandla  Base 

P.O.  Box  5800,  Albuquerque,  New  Mexico 

Attnt  Classified  DocumentOi vision 


Sperry  Gyroscope  Company 

Great  Neck,  L.I.,  New  York 

Attn:  Florence  W.  Tumbull,  Engr. 


Lib. 


library 


RCA  Defense  Electronic  Products 
Camden  2,   New  Jersey 
Attn:  Clarence  A.  Gunther,   Chief  Daf . 
DEP  -  Bldg.  IS,  F^oor  7 


Stanford  Research  Institute 

Document  Center 

Menlo  Park,   California 

Attn:  Mary  Lou  Fields— Requisitions 

Sylvanla  Electric  Products,  Inc. 
100  First  Avenue 
Waltham  5U,  Mass. 

Attnt  Charles  A.  Thomhill,   Reports  Lib. 
Waltham  Laboratories  Library 

Systems  Laboratories  Corporation 
15016  Ventura  Boulevard 
Sherman  Oaks,   California 
Attn:  Dr.   W.C.  Hoffman 

Technical  Research  Group 

17  Union  Square  West 

New  York  3,   N.Y. 

Attn:  M.L.  Henderson,   Librarian 

A. 3.  Thomas,  Inc. 

l6l  Devonshire  Street 

Boston  10,  Mass. 

Attn:  A.S.  Thomas,   President 


Engr. 


Westinghouse  Electric  Corp. 
2519  Wilkens  Avenue 
Baltimore  3,  Maryland 
Attn:  Kent  M.  Hack,  Qngr.  Ub. 

Wheeler  Laboratories,  Inc. 
122  Cutter  Mill  Road 
Oreat  Neck,  New  York 
Attn:  Mr.  Harold  A.  Wheeler 

Zenith  Plastics  Company 
Box  91,  Oardena,  California 
Attn:  Mr.  3.S.  Oleesky 

Library,  Oeophyslcal  Institute  of  the 
University  of  Alaska 
College,  Alaska 

University  of  California 

Berkeley  1:,  California 

Attnt  S.  Silver,  Prof.  Engr.  Science 

Division  of  Elec.  Engr. 

Electronics  Research  Lab. 

Electronics  Research  Laboratory 
University  of  California 
332  Cory  Hall,  Berkeley  U,  Calif. 
Attn:  J.R.  Whinnery,  Prof. 

California  Institute  of  Technology 
Jet  Propulsion  Laboratory 
U800  Oak  Grove  Drive 
Pasadena,  California 
Attn:  I.E.  Newlan 

California  Institute  of  Technology 
1201  E.  California  Street 
Pasadena,  California 
Attn:  Dr.  Charles  H.  Papas 

Carnegie  Institute  of  Technology 
Schenley  Park 

Pittsburgh  13,  Pennsylvania 
Attn:  Albert  E.  Heins 

Cornell  University 

School  of  Electrical  Engineering 

Ithaca,  New  York 

Attn:  Prof.  G.C.  Dalman 

University  of  Florida 
College  of  Eingineering 
Gainesville,  Florida 
Attnt  Prof.  M.H.  Latour, 

Engineering  Sciences  Lib. 

Georgia  Institute  of  Technology 
Engineering  Experiment  Station 
Atlanta,  Georgia 
Attn:  Mrs.  J.H.  Crosland 

Georgia  Tech.  Library 

Harvard  University 

Technical  Reports  Collection 

Gordon  McKay  Library,  303A  Pierce  Hall 

Oxford  Street,  Cambridge  38,  Mass. 

Attn:  Mrs.  E.L.  Hufsohmidt,  Lib. 

Harvard  College  Observatory 
60  Garden  Street 
Cambridge,  Mass. 
Attn:  Dr.  Fred  L.  Whipple 

University  of  Illinois 
Serials  Dept.  -  220S  Library 
Urbana,  Illinois 

Electrical  Engineering  Res.  Lab. 
University  of  Illinois 
Urbana,  Illinois 
Attn:  Anteima  Lab« 


The  Johns  Hopkins  University 
Homewood  Campus 
Depai^ment  of  Physics 
Baltimore  18,  Maryland 
Attn:  Professor  Donald  E.  Kerr 

The  Johns  Hopkins  University 
8621  Georgia  Avenue 
Silver  Spring,  Maryland" 
Attn:  Mr.  George  L.  Selelstad, 
Applied  Physics  Lab. 
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Massachusetts  Institute  of  Technology 
Research  Laboratory  of  Electronics 
bocument  Room  20B-221 
Cambridge  39,  Mass. 

Massachusetts  Institute  of  Tech. 

Lincoln  Laboratory 

P.O.  Box  73 

Lexington  73,  Mass. 

Attn:  Henry  Straus,  Document  Rra.  A-229 

Unlverilty  of  Michigan 
Electronic  Defense  Group 
Engineering  Research  Institute 
Ann  Arbor,  Michigan 
Attn:  J.A.  Boyd,  Supervisor 

The  Oniversity  of  Michigan 
Willow  Run  Laboratories 
WUlow  Rub  Alrpwt 
Ypsilanti,  Michigan 
Attn:  K.M.  Siegel,  Head 

Theory  and  Analysis  Dept. 

University  of  Michigan 
Willow  Run  Laboratories 
Engineering  Research  Institute 
Willow  Run  Airport 
Ypsilanti,  Michigan 
Attn;  Librarian 

The  University  of  Minnesota 

Minneapolis  lli,  Minnesota 

Attn:  Mr.  Robert  H.  Stumm,  Library 

Northwestern  University 
Microwave  Laboratories 
Evanston,  Illinois 
Attn:  Professor  R.E,  Beam 

Ohio  State  University 
Research  Foundation 
Columbus  10,  Ohio 
Attn:  Dr.  T.E.  Tice 

Dept.  of  Electrical  Engr. 

The  University  of  Oklahoma 

Norman,  Oklahoma 

Attn:  Prof.  C.L.  Farrar 

Chairman  Electrical  Engr. 

Polytechnic  Institute  of  Brooklyn 
Microwave  Research  Institute 
5?  Johnson  Street 
Brooklyn  1,  New  York 
Attn:  Vj*.  A.E,  Laemmel 

Syracuse  University 
Research  Institute 
Collendale  D-6 
Syracuse  10,  New  York 
Attn:  Dr.  C.S.  Grove, Jr. 

Director  of  Engr,  Research 

The  University  of  Texas 
Electrical  Engr.  Research  Lab. 
P.O.  Box  8026,  University  Station 
Austin  12,  Texas 
Attn:  J.R.  Oerhardt,  Assistant  Dir. 

The  University  of  Texas 

Defense  Research  Lab. 

Austin  12,  Texas 

Attn:  Claude  W.  Horton,  Phys.  Lib. 

The  University  of  Toronto 
Department  of  Electrical  Engineering 
Toronto,  Canada 
Attn:  Prof.  G.  Sinclair 

LTI  Research  Foundation 

P.O.  Box  709 

Lowell,  Mass. 

Attn:  Dr.  Charles  R.  Mingina 

University  of  Washington 

Dept.  of  Electrical  Engineering 

Seattle  5,  Washington 

Attnt  G,  Held,  Associate  Professor 

Polytechnic  Institute  of  Brooklyn 
Microwave  Research  Institute 
55  Johnson  Street 
Brooklyn  1,  Mew  York 
Attn:  Dr.  A. A.  Oliner 


Electronics  Division 
Rand  Corporation 
1700  Main  Street 
Santa  Monica,  California 
Attn:  Dr.  Robert  Kalaba 

National  Bureau  of  Standards 

Washington,  D.C. 

Attn:  Dr.  V.K.   Saunders 

Applied  Mathematics  and  Statistics  Lab. 
Stanford  University 
Stanford,  California 
Attn:  Dr.  Albert  H.  Bowker 

Dept.  of  Physics  and  Astronomy 
Michigan  State  University 
East  Lansing,  Michigan 
Attn:  Dr.  A.  Leltner 

Ion«phere  Research  Laboratory 
Pennsylvania  State  College 
State  College,  Pennsylvania 
Attn:  Prof;  A.H.  Waynick,  Director 

Institute  of  Mathematical  Sciences 

25  Waverly  Place 

New  York  3,  New  York 

Attn:  Mrs.  JoAn  Segal,  Librarian 

School  of  Electrical  Engineering 

Purdue  University 

Lafayette,  Indiana 

Attn:  Professor  F.V.  Schultz 

University  of  Tennessee 
Khoxville,  Tennessee 
Attn:  Dr.  Fred  A.  Ficken 

California  Institute  of  Tech. 
1201  E.  California  Street 
Pasadena,  California 
Attn:  Dr.  A.  Erdelyi 

VJayne  University 

Detroit,  Michigan 

Attn:  Professor  A.F.  Stevenson 

Mathematics  Department 
Stanford  University 
Stanford,  California 
Attn:  Dr.  Harold  Levine 

University  of  Minnesota 
Minneapolis  lli,  Minnesota 
Attn:  Prof.  Paul  C.  Rosenbloom 

Department  of  Mathematics 
Stanford  University 
Stanford,  California 
Attn:  Professor  Bernard  Epstein 

Applied  Physics  Laboratory 
The  Johns  Hopkins  University 
8621  Georgia  Avenue 
Silver  Spring,  Maryland 
Attn:  Dr.  B.S.  Oourary 

(2)Exchange  and  Gift  Division 
The  Library  of  Congress 
Washington  25,  D.C. 

Electrical  Engineering  Dept. 
Massachusetts  Institute  of  Technolnqy 
Cambridge  39,  Massachusetts 
Attn:  Dr.  L.J,  Shu 

Nuclear  DevelopiBent  Associates,  Inc. 

5  New  Street 

White  Plains,  New  York 

Attn:  Library 

Lebanon  Valley  College 
Annville,  Pennsylvania 
Attn:  Prof.  B.H.  Bissinger 

Dept.  of  Physics 
University  of  Pittsburgh 
Thaw  Hall 
Pittsburgh  13,  Pa. 
Attn:  Dr.  Edward  Gerjuoy 

Dept.  of  Physics 
Amherst  College 
Amherst,  Massachusetts 
Attn:  Dr.  Arnold  Arons 


California  Institute  of  Tech. 
Electrical  Engineering 
Pasadena,  California 
Attai  Dr.  Zofarsb  A.  JCaprielian 

Dr.  Rodman  Doll 
311  W.  Cross  St. 
Ypsilanti,  Michigan 

California  Institute  of  Technology 
Pasadena  b,  California 
Attn:  Mr.  Calvin  Wilcox 

Mr.  Robert  Brockhurst 

Woods  Hole  Oceanographic  Institute 

Woods  Hole,  Massachusetts 

National  Bureau  of  Standards 
Boulder,  Colorado 
Attn:  Dr.  R.  Gallet 

Mrs.  Jane  Scanlon 
22lt-03  67th  Avenue 
Bayside  dx,   New  York 

Dr.  Solomon  L.  Schwebel 

3689  Louis  Road 

Palo  Alto,  California 

University  of  Minnesota 
The  University  Library 
Minneapolis  lli,  Minnesota 
Attn:  Exchange  Division 

Professor  Bernard  lYiedman 
Dept.  of  Mathematics 
University  of  California 
Berkeley,  California 

Lincoln  Laboratory 

^ssachusetts  Institute  of  Tech. 

P.O.  Box  73 

Lexington  73,  Massachusetts 

Attn:  Dr.  Shou  Chin  Wang,  Rm.  C-351 

Melpar ,  Inc . 

3000  Arlington  Boulevard 

Falls  '^hurch,  Virginia 

Attn:  Mr.  K.S.  Kelleher,  Section  Head 

Hq.  Air  Force  Cambridge  Res.  Center 
Laurence  G.  Hanscom  Field 
Bedford,  Massachusetts 
Attn:  Lt.  Mark  J.  Beran 

Electronics  Research  Directorate 
Hq.  Air  Force  Cambridge  Res.  Center 
Laurence  G.  Hanscom  Field 
Bedford,  Massachusetts 
Attn:  Dr.  Philip  Newman,  CRRK 

Hq.  Air  Force  Cambridge  Res.  Center 
Laurence  G.  Hanscom  Field 
Bedford,  Massachusetts 
Attn:  Mr,  Francis  J,  Zucker 

Mr.  N.C.  Oerson 

Trapelo  Road 

S,  Lincoln,  Massachusetts 

Systems  Development  Corp» 
2L00  Colorado 
Santa  Monica,  California 
Attn:  Dr.  Richard  B.  Barrar 

Columbia  University 
Hudson  Laboratories 
P.O.  Box  239 
ll»S  Palisade  Street 
Dobbs  Ferry,  New  York 
Attn:  Dr.  M.W.  Johnson 

Institute  of  Fluid  Dynamics 

and  Applied  Math,  Univ.  of  Maryland 

College  Park,  Maryland 

Attn:  Dr.  Elliott  Montroll 

Dept,  of  Electrical  Engineering 
Washington  University 
Saint  Louis  5,  Missouri 
Attn:  Prof.  J.  Van  Bladel 
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Dept.  of  the  Navy 

Office  of  Maval  Res.  Branch  Office 

1030  Eo  Green  Street 

Pasadena  1,  California 

Braodeis  University 
Waltham,  Massachusetts 
Attn:  Library 

General  Electric  Comparer 
Microwave  Laboratory 
Electronics  Division 
Stanford  Industrial  Park 
Palo  Alto,  California 
Attn:  Library 

Smyth  Research  Associates 
3930  Uth  Avenue 
San  Diego  3,  California 
Attn:  Dr.  John  B,  Smyth 

Electrical  Engineering 
California  Inst,  of  Technology 
Pasadena,  California 
Attn:  Dr.  Georges  Weill 

Naval  Research  Laboratory 
Washington  2^,   D,C. 
Attn:  Dr.  Henry  J.  Passerini 
Code  ?278  A 

Dr.  George  Kear 
S  Culver  Court 
Orinda,  Calif. 

Brooklyn  Polytechnic 

85  Livingston  Street 

Brooklyn,  New  York 

Attn:  Prof.  Nathan  Marcuvits 

Dept.  of  Electrical  Engineering 
Brooklyn  Polytechnic 
85  Livingston  Street 
Brooklyn,  New  York 
Attn:  ^h,   Jerry  Shmoys 

Dept*  of  Mathematics 
University  of  New  Mexico 
Albuquerque,  New  Mexico 
Attn:  Eh*.  I.  Kolodner 
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